FRAME CONSTANTS OF GABOR FRAMES NEAR THE 

CRITICAL DENSITY 



A. BORICHEV, K. GROCHENIG, AND YU. LYUBARSKII 

Abstract. We consider Gabor frames generated by a Gaussian function and de- 
scribe the behavior of the frame constants as the density of the lattice approaches 
the critical value. 



1. Introduction 

In this article we study the stabihty problem for the expansions of functions 
on the real line with respect to a discrete set of phase-space shifts of a Gaussian, 
precisely 

f(x) = J2 Cfc^e^'^^'^'e-'^^"-''^)' . (1.1) 

k,l€Z 

Expansions of such form (with a = 1, 6 = 1) were introduced by D. Gabor in his 
classical article [5]. Now expansions of type (11. ip . so-called Gabor expansions, ap- 
pear in signal processing, quantum mechanics, time-frequency analysis, the theory 
of pseudo differential operators, and other applications. 

During the last decades an extensive theory of expansions (11.11) as well as more 
general Gabor expansions has been developed (see, for instance [3,6] and the ref- 
erences therein). However, not much is known about numerical stability property 
of such expansions. 

In modern language, the existence of Gabor expansions is derived from frame 
theory. To fix terminology and notation, take some g G L^(R), it will be called a 
window function, and let A = MZ^ C be a lattice in M^, where M is a 2 x 2 
invertible real- valued matrix. Given a point A = (x, ^) in phase-space M^, the 
corresponding time-frequency shift is 

The set of functions Q{g,A) = {n\g : A G A} is called the Gabor system generated 
by g and A. We say that such a system is a Gabor frame or Weyl-Heisenberg frame, 
whenever there exist constants A,B>0 such that, for all / G L^(]R), 

MfWlHR) < E K/'^A^?)l^(R)P < . (1.2) 

AeA 
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The (best possible) constants A = A{A,g) and B = B{A,g) in fll.2p are called the 
lower and upper frame bounds for the frame Q{g,A). 

There is a standard procedure for constructing expansions of type (11 .ip for each 
Gabor frame. Namely there exists a dual window 7 G L^(R), such that every 
/ G L^(]R) can be expanded as a Gabor series 

f = J2(f,^M)^x9- (1-3) 
aga 

Such dual window is, in general, non-unique. We refer the reader to e.g. [6] for an 
exposition of Gabor analysis and related matters. 

The property of the system Q{g,A) to form a frame in L^(M) depends (among 
other factors) on geometrical characteristics of A. We say that the area of its 
fundamental domain s{A) = Area(M[0, 1)^) = |detM| is the size of A. By the 
density of A we mean d{A) = s(A)~^; for the lattice case this definition of density 
coincides with numerous standard density definitions (see e.g. [16]). 

We refer to [9] for a comprehensive account of the density theorems for Gabor 
frames. In fact for any window function g the condition s{A) < 1 is necessary for 
Q{g,A) to be a frame in L^(R). For "nice" windows g (in the Schwartz class, say) 
a fascinating form of the uncertainty principle, the so-called Balian-Low Theorem 
(BLT), requires even that s(A) < 1 for ^(^f. A) to be a frame [2]. 

The results of [13,17] yield in particular that in the case of the Gaussian window 

go{t) teR. 
the condition s(A) < 1 is also sufficient: 

Theorem A. The set Q{go,A) is a frame for L^(]R) if and only if s{A) < 1 . 

Together with BLT this implies that the lower frame bound A = A{A) must 
tend to 0, as the size of the lattice s(A) approaches one. Thus the original Gabor 
series (II. ip with a = 1, 6 = 1 corresponds to the critical case s{A) = 1 and 
does not provide an L^-stable expansion. In this case, there exist L^-functions 
with polynomially growing coefficients. See [10, 15] for the convergence properties 

of dni). 

In this article we are concerned exclusively with Gabor frames for the Gaussian 
window Q{go,A) for the square lattice A(a) = aZ x aZ and study the behavior of 
its frame constants A{a) = A{A{a)) and B{a) = B{A{a)) near the critical density 
d{A) = 1. The main result of the article reads as follows. 

Theorem 1.1. There exist constants < c < C < 00 such that for each a G 
(1/2, 1) the frame bounds A{a), B{a) for the frame Q{gQ,A{a)) satisfy 

c{l - a^) < A{a) < C{1 - a^) (1.4) 

and 

c < B{a) < C. (1.5) 

""^In his article [5] Gabor considered expansions of functions / which possess additional decay in 
time and frequency. As we now know (see e.g. [15]), for such functions the series (|l.ip converges 
in L^{M.). 
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Remark 1.2. A similar statement holds for arbitrary rectangular lattices. The 
values of c, C in this theorem then depend upon the shape of the lattice. Neverthe- 
less, one can prove that given a number i^' > there exist constants c and C valid 
for all matrices M such that the diameter of the fundamental domain M(0, 1]^ does 
not exceed K. 

The ratio B{K) / A{h) plays the role of the condition number for the frame 
Q{gQ,h). Thus Theorem 11.11 says how fast does the frame Q{gQ,h) "numerically 
degenerate" as its density approaches the critical value. 

The asymptotical behavior A{a) x {1 — a?) has been first observed numerically 
by Thomas Strohmer [19] and by Peter Sondergaard [18]. Moreover, the numerical 
simulation in [19] allowed us to guess the construction which gives the second 
inequality in fll.4p . This construction is described in Section 4 below. 

Next let gi{t) = (cosh7r7t)~^, 7 > 0, be the hyperbolic cosine function. Janssen 
and Strohmer [12] have shown that Q{gi, aL x is a frame, if and only if ah < 1. 
To do this, they showed that the frame bounds for Q{gi,aZ x 6Z) are equivalent 
to those of Q{go, o-^ x bZ) with the Gaussian go and applied Theorem A. Therefore 
we obtain the same asymptotic estimates for the frame bounds for the hyperbolic 
cosine. 

Corollary 1.3. There exist constants < c < C < 00 such that for each a G 
(1/2, 1) the frame bounds A{a), B{a) for the frame Q{gi, A{a)) satisfy 

c{l-a^) <A{a) <C{l-a^) 

and 

c < B{a) < C. 

The proof of Theorem 11.11 involves both time-frequency methods and methods of 
complex analysis. We use complex analysis in order to obtain the upper estimate 
for A{a) and the Gabor analysis in order to obtain the rest of the statements in 
Theorem 1.1 (though a pure complex-analytic proof is also available). In particular 
we apply Walnut's estimates for the norm of the frame operator [20], and also 
precise decay estimates for the dual window established in [7]. The upper bound 
A{a) < C(l — a^) will be established by the construction of a concrete example. 
We produce a function fa (depending on the lattice A(a)), such that 

E \{fa,^x9o)\'<C{l-a')\\fJl,^^y 

AeA(a) 

By using the Bargmann transform, we translate our problem into one of finding 
entire functions in the Bargmann- Fock space whose restrictions to A(a) are "small" 
with respect to their Fock norms. 

The paper is organized as follows. In the next section we discuss the estimates 
for B{a). Furthermore, we give the lower estimates for A{a). Here we mainly follow 
the arguments from [7]. In section 3 we recall the definition of the Fock space J-' of 
entire functions, and discuss the relations between the frame property of the system 
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Q{go, A (a)) and sampling in J^. We also recall basic properties of the Weierstrass a- 
function. In section 4 we use these facts and also special " atomization" techniques 
in order to construct the example which delivers the upper estimate in fll.4p . 

Notation: To avoid dealing with too many intermediate constants, we use the 
standard notation f -< g to express an inequality f{x) < Cg{x) for all x with a 
constant C independent of x (and possibly other parameters). Likewise, f ^ g 
means that there exist A, i? > such that Af{x) < g{x) < Bf{x) for all x. 

2. Time-Frequency Methods To Estimate Frame Bounds 

The estimate (11. 5p on the upper frame bound B{a) can be obtained in various 
ways. In particular, we can use Walnut's estimates, which give a sufficient condition 
for the Gabor frame operator to be bounded [20]. This result also follows from the 
Polya-Plancherel type inequalities for functions in the Bargmann-Fock space, see 
below Section 3 for more details. 

To obtain the lower estimates for A{a) we need to show the invertibility of the 
Gabor frame operator and to estimate the norm of the inverse operator. We will 
approach this problem by using information about a suitable dual window 7 and 
then apply Walnut's estimates to the Gabor expansion (11.31) . 

To state Walnut's result we need the following definitions. 

Let W be the the Wiener amalgam space of functions on the real line defined by 
the norm 

hWw = sup \g{t + k)\ . 
fcez *e[o,i] 

Given a function g in L^(]R), consider the Gabor system Q{g,A) and the corre- 
sponding synthesis operator -Dc,,a, 

Dg^AC = "^cxnxg, 

AeA 

and the analysis operator Cg,A, 

iCg,Af)iX) = {f,nxg), AG A. 

If Dg A acts continuously from £^(A) to L^(]R), then Cg^A acts continuously from 
L\R) to i\A), and Cg^A = ^*,a- 

The following lemma from [20] gives an estimate for ||Z}g^A||/2-^L2: 

Lemma 2.1. If g G W and A = aZ^, then Dg A is bounded from £^(A) to L^(]R) 
and 

||^g,A||L2(i;) < {1 + a''^)\\g\\w 

Since, obviously, in our situation B{a) = ||Cg,A|li2^;2 = 11-0^^11^2^2,2, we obtain 

Corollary 2.2. If g e W and a> 0, then 

B{a)<{l + a-'Y\\g\\l.. (2.1) 

To treat Gabor frames with Gaussian window, we need to evaluate the amalgam 
space norm of functions with Gaussian decay. 



FRAME CONSTANTS NEAR THE CRITICAL DENSITY 5 

Lemma 2.3. Assume that K>0, |7(t)| < e"^'^* . Then 

hllvK < 2 + (2.2) 
Proof. For n > 1, n G Z, we have 



sup |7(n + t)| < e"™' < / e""''*' rft 
e[o,il A-i 



te[o,i] 

and likewise for n < — 1, n G Z, we have 

/•|n|-l 

sup |7(n + t)| < e-"''(l'^l-i)' < / e~"'^*'dt. 

te[0,l] J|n,|-2 

Consequently, 



V sup \-^{t + n)\ <2+ e~ 
„^r, te[o,i] Jm 



TTKt' 



dt = 2 + K~^/^ . 



□ 



As a consequence we obtain an estimate on the upper frame bound of Gaussian 
Gabor frames. 

Proposition 2.4. The upper frame bound B{a) of Q{go.,a7?), 1/2 < a <1, satis- 
fies the estimate 

1 < B{a) < 100. 

Proof. For the upper estimate we use (12. ip and (12.21) with k = 1. 

To get the lower estimate we consider the sum (11.21) foT f = g = qq. Then 

|(^o,7rA^o)P > \\9of, 

AeA(a) 

which yields the desired estimate. □ 

The time- frequency methods also yield the lower estimate in (11.41) . This estimate 
requires the existence and some knowledge about a dual window. If Q{g,A) is a 
frame, then by the frame theory there exists a dual window 7 G L^(]R), such that 
every / G L^(M) possesses a(n unconditionally convergent) series expansion (Gabor 
expansion) of the form 

/ = 7^X9)7^X1 = D^^ACg^Af . 

AgA 

For the square lattice A(a), Lemma [2.11 yields the following bound: 

\\f\\hm<{a-' + m\irwJ2\^f,7T,g)\\ 

AeA 

Consequently, the lower frame bound 74(a) can be estimated from below as 

A(a)>((a"i + l)2||7||^)-'. (2.3) 
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Proposition 2.5. For the square lattice A(a), 1/2 < a < 1, the lower bound A{a) 
of the Gaussian frame Q{gQ, A{a)) obeys the estimate 

A{a) y l-a"^. 

Proof. In [7] the authors consider the Gaussian Gabor frame Q{gQ, A{a)). For this 
frame they construct a dual window 7 such that 

|7(t)| <Ce-"''*' 

with K X 1 — a^. 

By Lemma [2.31 we have 

hik ^ 2 + ^ (1 - a2)-i/2, 

and the desired estimate follows now from (12.31) . □ 

3. Complex Methods 

3.1. Fock space. We recall the definition and basic properties of the Fock space. 
We refer the reader to [4], [6] for detailed proofs and also for a discussion of nu- 
merous applications of this space to signal analysis and quantum mechanics. 
The Fock space T is the Hilbert space of all entire functions such that 

\\F\\:= [ |F(2)|V^I"l'rfm, < 00, 
Jc 

where druz is Lebesgue measure on C. 

The natural inner product in J-' is denoted by (■, 
We will use the following well-known facts: 

(a) The point evaluation is a bounded linear functional in JF, and the corre- 
sponding reproducing kernel is the function w 1— > e'^^'^, i.e., 

F(z) = (F,e-^"-)^, FG^. (3.1) 

(b) One defines the Bargmann transform of a function / G L^(]R) by 

/ ^ Bf{z) = F{z) = [ f{t)e-^'\^^''dt. 

Jr 

The Bargmann transform is a unitary mapping from L^(]R) onto JF. 

(c) In what follows we identify C and M^. In particular for each ( = + irj E C 
we write tt,^ = 7f{^,r)) ■ Define the Fock space shift : JF ^ JF by 

f3^F{z) = e^"«''e-"l^l'/2e"^"F(;z - C). 

Then is unitary on JF, and the Bargmann transform intertwines the Fock space 
shift and the time- frequency shift: 

(3^B = Bn^. (3.2) 

(d) 

Bgo = 2-l/^ (3.3) 
here as above go is the Gaussian function. 
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(e) It follows from and that 

Taking into account the reproducing property (13. ip . we can rewrite the frame 
property (11.21) of Q{go,A) as the sampling inequality 

AeA 

In the case of square lattice, A is symmetric with respect to the real line, and we 
have 

v4||F||^ < -^|F(A)|V"I^I' < FeJ^. (3.4) 

AeA 

(f) Let 1/2 < a < 2, and let w E C, w 0. Consider the entire function 
^aA^) = e"'^"'/"'. Then 

|$a,«;(^)| X e'^l^l', \Z-W\<1. 

This statement can be checked by direct inspection. 

3.2. Reformulation of the main result. The remaining part of Theorem 11.11 
can now be reformulated as follows. 

Theorem 3.1. Let A(a) = aZ^, and let A{a) be the best possible constant in the 
left hand side inequality of (13.41) . Then for 1/2 < a < 1 we have 

A{a) -< 1 - 

To prove this theorem we need to find a constant K and functions F = Fa E 
such that 

K{l-a^)\\Faf> J2 l^a(A)|V-l^l'. (3.5) 

AeA(a) 

3.3. Weierstrass a-function. The construction of the functions Fa in the next 
section is motivated by the properties of the classical Weierstrass cr-function. Let us 
recall its definition and basic properties. We refer the reader to [1] for a systematic 
study of this function and also to [7] for its applications in Gabor analysis. 

Given a lattice A C C we denote 

a{A,z)=z n (l-f)^*^^^^' 

AeA\{0} 

This product converges uniformly on compact sets in C to an entire function with 
A as the zero set. This is a function of order 2; moreover there exists c/a G C such 
that 

|a(A, z)e'^''''\ X ei^(^)"'l^l', dist(A, z) > e > 0. 

Here s(A) is the area of the fundamental domain of A. See [8] and also [7]. 
Once again, let A (a) = aZ^. A direct inspection shows that d^a) = 0, so that 

|a(A(a),^)| xe^"~'l"l', dist(A„z) > £ > 0. (3.6) 



V 2 
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This relation allows one to mimic the weight function e'^'^'^''^ in the definition of 
the Fock space by the absolute value of an analytic function. 

4. Proof of (1331) 

4.1. Explicit construction. If a is in a compact subinterval of (1/2, 1), one can 
take Fa = 1 and obtain fl3.5p with some appropriate constant K. Therefore, from 
now on, we assume that a is sufficiently close to 1, say 0.999 < a < 1. Given such 
a, we take R = R{a) such that 

2(1 - a") < < _ 

and 

riR := 7r(l - R-'^'^)R^ G N. 
We need some additional notation: 

h^ = l-R-^l\ 
Cm,n = b~^{m + in), 

Qm,n = {x + iy e C:\x- b'^m\ < \y - b'^n\ < b-^/2}, 

Dn = {z eC:\z\< R}, 

D'r = U{Qm,n ■■ \Cm,n\ < R ' 3} , D'^, = Dr \ D'j,, 

Afn = {{m,n) eZ^ -.Q^^^cD'j,}, 
qr = Card Mr, pr = ur - qr . 

We have 

{z : R-l <\z\ <R} G D'^G {z : R- 4:<\z\ < R}. 
Using the appropriate segments of radii of the disc Dr we split D'^ into the "sectors" 

PR 

D'k=[J Ak, 

k=l 

such that 

m < diam Ak < M, Area Ak = 
for some m, M independent of a. Denote the center of mass of A^ by 

Ck = b^ [ Cdm^. (4.1) 

We can find c independent of a such that 

{w: \w - Ck\ < c} C Ak, k = l,...,pR. 
We are going to verify the estimate (13.51) for the function 

\ PR 

z 



(m,n)eA/'fl\(0,0) 
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The zero set of the function Fa is 

Za = {Cm,n : |Cm,n| < i? - 3} U {Ck]Hi • (4.3) 

By construction, the total number of zeros of Fa is = nR^lP'. 
In order to prove (13. 5p . we need to estimate both ||-Fa||3r and 

ll^a|la:= E |i^a(a(m + zr2))|V™'("^'+"'). 

4.2. Estimate of ||-Fa||^- To estimate the norm of Fa in the Fock space, we com- 
pare the logarithm of the modulus of the polynomial Fa to a subharmonic function 
ur whose growth is easy to control. 
Consider the subharmonic function 

£^^. = 1 iLf; , , . N<^' 

An easy estimate shows that 



ur{z) = / log 
'\C\<R 



1 log 1^1 -7ri?2 log + \z\>R. 



ur{z) < —\z\'^, \z\ > R. 



2 

We use the following approximation lemma. 
Lemma 4.1. For each e > there exist constants < c^e) < C{e) < oo such that 
c{e)\Fa{z)\ < e'''"^^^) < C{e)\Fa{z)l dist(^,Zj > e. (4.4) 

and 

c{e)\Fa{z)\ < e'"^^^'\ dist(^,Zj < e. (4.5) 

Remark. Since the set M is invariant with respect to rotation by 7r/2 around the 
origin, we find that 

r — = F ' ^0. 

^-^ m + in ^ (m + inV 

So the first factor on the right-hand side of (14.21) in the definition of Fa can be 
written as 

z 



c 

(m,n)6A/'H\(0,0) 

{m,n)eNR\{0,0) 

Consequently, the function can be viewed as a truncated version of the Weier- 
strass (T-function and estimates fl4.4l) and fl4.5p correspond to the growth estimate 
(13. 6p for the Weierstrass a-function. 

We postpone the proof of this technical lemma until subsection 14.41 Assuming 
that Lemma [4.11 is already proved, an estimate of ||-Fa||3^ is straightforward. 



Lemma 4.2. 



\Fa\\l y R'/' - {1 - a^)-\ 
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Proof. Let fl = {z E C : \z\ < R, dist{z, Za) > 1/10}. Using Lemma WTi we find 
that 



Jn 



\z\ = r 



We use that l — b'^ = R Furthermore, for every 1 < r < i?, the circle z : 
intersects with Q on at most half of its length. Therefore, 

fR , p3/2 /-Ri/^ 

J(a, R) y / e-^'''"'"'"tdt = 1— / e-^^du x R'/\ 

and the statement of the Lemma now follows. □ 
Remark. A similar argument shows that 

|F,(^)|V-l^l'dm,^0, (4.7) 



\z\>R-4 

as a — > 1 or equivalently, as R ^ oo. 

4.3. Estimate of \\Fa\\l. 

Lemma 4.3. For Fa as in f l4.2p we have 

\\Fa\\l = J2 \Fa{a{m + zn))|2e-™'("'+"') x L (4.8) 

m,n 

Proof. We have 

ll^-lla=( E + E )|i^a(a(m + 2n))|V™'('"'+"') = Si(a) + S2(a). 

(m,ri)GAfl (m,ri)0Afl 

In order to estimate Si (a), we observe that Fa{Cm,n) = 0, for (m, n) G A/}? and 
|Cm,n ~ Cii"^ + ^'^)| = {b^^ — a)\m + in\ < 2R~^^'^\m + in\. 
Since |m + m| < i? for (m, n) G A/r, and a is sufficiently close to 1, we have 

\Cm,n- a{m + in)\ <l, {m,n) e Ur. 

5 

Denote -Dm,n = G C : \z — Cm,n\ < 1/4}- By part (f) of subsection 13.11 there 
exists a function $m.n(^) that is holomorphic on Dm,n and satisfies 

|$m,n(^)|xe'''fl^^ ^GD^,„. (4.9) 
Then for each (m, n) G A/}? the function 

iz) = ^^ 

is holomorphic in ii'm.n and possesses the properties 

^m,n{Cm,n) = and |^ 

m,n 

By Cauchy's theorem, the functions ^1/^,^ uniformly bounded on -D^_^ = 
{z G C : I2; — Cm,n| < 1/5}, and hence 

|^„^„(a(m + m))| ^ |(a - 6~^)(m + m)| ^ R~^''^{m^ + n^)^/^ (m,n) G A/^j. 
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Returning to the function Fa and using (14 .Qp once again, we obtain 
The mean value inequahty for \(^rn,n{0'ij^ + in)Y^^^^ ^^1 now yields 



Since all discs -Dm,n are disjoint we obtain 

J\z\<R Jo 

Finally, for arbitrary {m,n), the mean value theorem yields 

and, by (14.71) we obtain 

S2(a) ^ / |F„(z)pe-''l^l'(im^ ^ 0, as a / 1. (4.11) 

J\z\>R-4: 

The estimates (l410|) and (I4TTD yield 

= 5^|F„(a(m + m))|V™'("'+"') -< 1. (4.12) 

ni,n 

The opposite relation follows from Lemma Upland from the lower estimate on A{a) 
estabhshed in Proposition 12.51 □ 

Relation (13. 5p follows immediately from Lemma 14.21 and the estimate (14.80 (or 
even ( Km ). 

4.4. Proof of the approximation lemma. The proof of Lemma [4. II is based on 
atomization techniques, see e.g. [14]. First we rewrite (14.41) and (14.51) in an additive 
form. We must prove that 

\og\Faiz)\=b\Riz) + 0{l), dist(z,Zj >£, (4.13) 

and 

log|F,(2;)| < b\R{z) + 0{l), dist(2;,Z„) < e, 

where Za = {Cm,n '■ {m,n) G A/r} U {Ck}k=i &s in (14.30 . and the quantities 0(1) in 
the right-hand sides of these relations are bounded uniformly with respect to all 
a e (0.999, 1) and depend only on e. It suffices to prove (I4.13p . the second relation 
will then follow by the maximum principle applied to Fa^~l^ where ^a,w is defined 
in part (f) of subsection 13.11 
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Let Vr be defined by M . 



vr{,z) =b^ log 



1 - 



c 



dm, 



and let 



We have 





Z 


/ log 









PR / \ 

dm^, Wn{z) = n (^1 - I j • 



\og\Faiz)\ - b\R{z) = 

(log \Vr{z)\ - vr{z)) + (log \Wr{z)\ - wr{z)) = &^{R, z) + &2{R, z), (4.14) 

and we estimate separately each summand in the right-hand side of fl4.14p . 
Let (iisi{z,Za) > e. We have 



&^{R,z) = \og\VR{z)\-VR{z) 



Qo,o 



{m,n)£j\fR\{0,0}, (m,n)eJ^R\{0,0}, 

dist(2 •,Qm,,n )<io dist(2 )>io 



log \z\ — lo, 
log 



dlTLc 



1 



z 



Cm,r, 



log 



dm(^ . 



Jm,n 



It suffices to estimate just the second sum in the right-hand side because the 
first sum contains only a finite number (at most 1000, say) of uniformly bounded 
terms, and the first integral is bounded uniformly in z. 

Denote L{Q = log(l — z/Q. We then have 



Jm,n 



3? 



{HO - L{Cm,n)dm(;) 



We apply the second order Taylor expansion with the remainder term in the integral 
form: 



L(C)-L(C„,„) =L'(C„,„)(C-Cm,n) + ^i^"(Cm,n)(C-C™,n)' + ^ ^ L'" (s) {( - sfds . 



and use the fact that 



Qm.,Ti 



(C - Crn,n)dmc_ 



(C - C,m,nfdmc_ = 0. 



Then 

\jm,n I 



Qm^n 



1 



[s — zY 



— ds dm, 



-< 



+ 



dist{z,Qm,nY dist(0,Q 



which implies that 



G^{R,z) = 0{l). 
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Finally, 



\WR(z)-wn(z) 



dist(z,yis.)<M+io dist(2,yifc)>Af+io 



log 



z 



log 



with M as in (14. ip . The first term in the right-hand side contains just a finite 
number of summands and is always bounded. In order to estimate each ik from 
the second term we use the Taylor formula (now of the first order) with the same 
function L{() = \og{l—z/(). The choice of (k in (14. ip implies that J^^{(—(k)dm(; = 
0. Arguing as above, we obtain 

1 1 

^ dist(z, A)2 + dist(0,Afc)2' 

whence 

&2{R, z) = 0{l). 
This completes the proof of Lemma I4.1[ □ 
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